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Recent experimental results indicate that the dark matter sector may have a non-minimal struc- 
ture with a spectrum of states and interactions. Inelastic scattering has received particular attention 
in light of DAMA's annual modulation signal. Composite inelastic dark matter (CiDM) provides a 
dynamical origin for the mass splittings in inelastic dark matter models. We show that higher di- 
mensional operators in the CiDM Lagrangian lead to an admixture of inelastic and elastic scattering 
in the presence of parity violation. This scenario is consistent with direct detection experiments, 
even when parity violation is nearly maximal. We present an effective field theory description of 
such models and discuss the constraints from direct detection experiments. The CiDM model with 
parity violation has non-trivial phenomenology because of the multiple scattering channels that are 
allowed. 



I. INTRODUCTION 

Recent direct and indirect searches for dark mat- 
ter hint that the dark sector may have non-minimal 
structure and interactions. This is in sharp contrast 
with the standard scenario of weakly interacting mas- 
sive particles in which the dark matter is the lightest 
neutral state in a spectrum and interacts elastically off 
of Standard Model (SM) particles. The results of the 
DAMA experiment provide an example of an anomaly 
that challenges the standard dark matter picture [TJ |5] . 
In particular, if DAMA's measured annual modulation 
arises from inelastic dark matter (iDM), then DAMA 
can be reconciled with all other null results from direct 
detection experiments [21 H]- The presence of multiple 
states that lead to inelastic interactions may indicate 
novel dynamics in the dark sector. For example, iDM 
requires an 0(100 keV) splitting between the dark mat- 
ter states, which may be a sign that dark matter is com- 
posite [siEiiTiiiiiiiniiiiiinjiniiii]. 

Composite inelastic dark matter (CiDM) is a recent 
proposal that provides a dynamical origin for the 100 
keV mass splitting 5 . CiDM models have a ground 
state degeneracy that is split by the hyperfine interac- 
tion. In [5 , a minimal CiDM model was proposed where 
a new strong gauge group confines at low energies and 
the quarks that are charged under the new strong gauge 
group form "dark hadrons" after confinement. The inter- 
actions between the Standard Model and the dark sector 
are mediated by a kinetically- mixed C/(l)d- 

The parity of the U(l)d current determines how vis- 
ible the dark matter is to the Standard Model. The 
minimal CiDM model conserved parity and the axial- 
vector current interaction led to only inelastic inter- 
actions. However, parity can be explicitly broken by 
the dynamics of the new strong gauge group through a 
term. Parity violation leads to elastic interactions 



that arise from charge-radius scattering and are phe- 
nomenologically different from typical elastic dark mat- 
ter interactions because the charge-radius scattering is 
suppressed at low nuclear recoil energies. The recoil 
spectrum of this model looks strikingly similar to stan- 
dard iDM models due to the suppression from low re- 
coil energy interactions. CiDM models provide a new 
framework for studying kinematic scenarios where sev- 
eral types of scattering events are allowed. 

This article presents the theory of CiDM models 
with parity violation. Sec. ITT] describes the low energy 



III 



effective theory in terms of "dark mesons." Sec 
computes the direct detection phenomenology using a 
global fit while marginalizing over the uncertainty in the 



dark matter velocity distribution function. Sec. IV dis 



cusses constraints arising from QED tests and summa- 
rizes prospects for collider searches. 



II. MODELS OF CIDM 

In this section, the effective field theory for CiDM 
models is reviewed and the consequences of parity vio- 
lation in the new strong sector is explored. The high 
energy theory is a two flavor SU(N C ) gauge theory with 
a Lagrangian given by 

£ = £sM + ^CiDM (1) 

£ C iDM - --TrG 2 dflv + y L ip* L + * H iP*H 

where a = L,H , are Dirac fermions that are funda- 
mentals under the strong gauge sector and is the 
SU(N C ) gauge field strength. In Atomic Inelastic Dark 
Matter (AiDM) [6], the strong gauge sector is replaced 
by an Abelian gauge group where ^ a are charge ±1, 



respectively. If have chiral gauge charges, then the 
m a arise through a symmetry breaking interaction (i.e., 
m a = V a {<!>))■ If the theory is asymptotically free as in 
[5], then the theory will confine at a scale Ad- The bound 
states will be approximately Coulombic if > and 
the resulting spectrum is qualitatively similar to AiDM. 

To have the appropriate relic density the dark mat- 
ter must either be very heavy ( > 30 TeV) or the dark 
matter density must be generated non-thermally, pos- 
sibly linked to baryogensis. O(30 TeV) inelastic dark 
matter is not compatible with CDMS's null results 
QH [EEJ [T7] if it fits the DAMA signal and therefore the 
relic abundance of the dark matter needs to be gener- 
ated non-thermally. Assuming that there is a cosmolog- 
ical asymmetry generated early in the Universe between 
heavy quarks and light anti-quarks, 



n H -n a = -n L +n L ^0, 



(2) 



the dominant component of the dark matter will be in 
dark mesons with a single heavy quark [5]. 

The dark matter is the ground state of a ^h^l 
meson and will be denoted as n^. The 7Td is a spin 0, 
complex scalar with parity —1. Dark matter scattering 
is primarily a transition to the complex, spin 1 meson, 
Pd- The parity of this state is 



Ppdn = (-1)% ("I)" = 




(3) 



The mass splitting between the ltd and the pd arises 
through the hyperfine interaction and is suppressed when 
mji ^> wif,, Ad- In particular, 



kAa 



5m = 



N c m H 



m L < A d 



m L > A c 



(4) 



where Ad = N c g^/ATT is the 't Hooft coupling with N c = 1 
applying to Abelian gauge groups and k is an 0(1/N C ) 
constant. Note that m Vd ~ mjj in the heavy quark mass 
limit. For mass splittings 0(100 keV) and a dark matter 
mass near the weak scale, the confinement scale is ~ 100 
MeV. 

The dark matter mesons interact through a massive 
spin 1 gauge field A^ that kinetically mixes with t/(l)y 
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where D^cf) = d^4> ~ ^WdA^cj). After (j> acquires a vev, 
fy, Ad becomes massive and the mixing between the 



dark sector and the Standard Model can be diagonalizcd. 
The e mixing between the hyperchargc field strength and 
Z7(l)d is the source of the interactions between the Stan- 
dard Model fermions and the dark matter. Assuming 
that 



m Ad = 2g A f 4> < m z o 



(G) 



after electroweak symmetry breaking, the couplings can 
be diagonalized. The interactions relevant for dark mat- 
ter scattering are given by 



Ant = (Jd + ec e^EM) 



(7) 



where cq = cos W , Jem is the electromagnetic current 
and Jd is the current of the dark quarks. A more com- 
plete analysis of the interactions is given in Sec. |IV| 
Anomaly cancellation restricts the charge assignments 
of the two dark quarks leaving only three anomaly-free 
possibilities for the current of the dark quark sector. 



A. Axially Charged Quarks 

The types of interactions that are allowed depend 
on whether Jd is an axial or vector current. In [5] and 
[B], Jd is an axial vector current. The only anomaly- free 
axial charge assignment in terms of Weyl spinors is 
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(8) 



where the Dirac spinors ^ a = (ip a ,ip^). The masses of 
the quarks arise from Z7(l)d breaking through the Higgs 
mechanism 



h.c. 



(9) 



Because both the mass of the quarks and the Ad arise 
from the vev of <f>, there is a hierarchy between the gauge 
couplings and the Yukawa couplings 



25d 
yn ' 



(10) 



Fitting DAMA requires m nd ~ 100 GeV, while m Ad can 
in principle take on a range of values from 10 MeV to 
100 GeV. Eq. [lO] implies that the gauge coupling for the 
axial sector may be small in comparison to the Standard 
Model gauge couplings because yn is capped by pertur- 
bativity at 0(1). 

The effective operators describing the interactions 
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of the ltd — pd system are 



C. Parity Violation 
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The operators with coefficients denoted by d are sup- 
pressed by a factor of the relative velocity v Ie \, while 
the ones denoted by c are not velocity suppressed. The 
elastic scattering operator for the ir d is dimension 6 and 
velocity suppressed, resulting in an overall suppression 
of the elastic to inelastic scattering rate of v^ el . 



B. Vectorially Charged Quarks 

There are two anomaly-free charge assignments for 
vectorially charged dark matter: one gives the composite 
dark matter a charge and the other leaves it neutral. 
Charged dark matter will have an enormous scattering 
rate and will look qualitatively similar to standard elastic 
dark matter. The charge assignment that leaves the dark 
matter neutral will only scatter off higher moments of 
the charge distribution and will be suppressed at low 
recoil energy. The charge assignments for the neutral 
dark matter theory are 



Neutral Vector 
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(12) 



With these charge assignments, A d couples to a vector 
current and the allowed operators are 
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Vector off 



+ -4(^3^ + p^d^p Av )d v F^ 
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(13) 



d denotes operators that are velocity suppressed and c 
denotes unsuppressed operators. The leading operator 
that is not velocity suppressed is the elastic charge-radius 
operator, but this is a dimension 6 operator. Recent 
work on form factor-suppressed inelastic transitions in- 
dicates that this type of scattering may be an explana- 
tion for DAMA [24] [25] . The primary difference between 
form factor elastic scattering dark matter and iDM is the 
existence of a threshold in iDM. The next section illus- 
trates that it is possible to have dark matter dominantly 
scatter inelastically and have a residual form factor elas- 
tic contribution. 



In the two models above, parity determined the in- 
teractions of the dark meson fields; however, parity is 
not a fundamental symmetry of nature. If parity is bro- 
ken, both charge-radius scattering and inelastic scatter- 
ing are allowed without a velocity suppression. This is 
quite natural in strongly coupled CiDM models because 
CP violation arises from the dynamics of the strong sec- 
tor through the term 



Lf = 6 d Tr G d Gd, 



(14) 



and results in mixing between states of different parity. 
The size of 0d is not necessarily related to the size of 
0qcd a- n d in principle 0d could be 0{\). 

Because the Od term is a total derivative, its ef- 
fects only appear non-perturbatively. The dominant 
effect of the CP violation is to cause a small mixing 
between states of different parity. In QCD, for exam- 
ple, the 7T° with I G (J P ) = 1"(0") and the a with 
I (J p ) = 1~(0 + ) mix in the presence of a 0qcd term. 
A similar process will happen in the dark sector. When 
m L ;$ Ad, the mixing angle between fields of opposite 
parity is given by 



sin 9f ~ Od 



A d ' 



(15) 



The mixing vanishes in the limit where mj, — > because 
the 0d term can be removed by a chiral rotation of the 
If mL 3* Ad, the mesons form Coulombic bound 
states and the mixing angle is given by 



sin a v> = 



e d A d 

\ 2 d m L 



(16) 



where the matrix elements of the perturbing CP- 
violating Hamiltonian is set by the non-perturbative 
scale where the effects of 0d are not exponentially sup- 
pressed. As oii — > co, the CP- violating effects decou- 
ple and parity violation vanishes. Therefore, even if 
Od ~ 0(1), its effects on the interactions of the dark 
mesons might be small if ttil — > 0,oo. Maximal parity 
violation occurs when rriL — Ad- 

With an axially coupled U(l)d, the 7Td — aod inter- 
action becomes an elastic charge-radius operator with 
parity violation: 



Ar d a 0d - TT 7r d < ^ a Od<9 I /-F'd ' 
yv d 
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Therefore, the effects of 6d can be estimated by replac- 



3 



ing the field strengths in Eq. 11 and Eq. 13 with 1 



cos 26 f F^ + sin 26 r F ( 



(17) 



Therefore, turning on parity violation in the strong sec- 
tor allows admixtures of vector and axial vector inter- 
actions. The ratio of elastic to inelastic cross sections 
becomes a free parameter. The next section will show 
that an upper bound of Of < 0.08 is necessary to avoid 
direct detection constraints assuming that all c, d ~ 0(1) 
and itin ~ m 7Td . 



III. DIRECT DETECTION PHENOMENOLOGY 

Novel features in the direct detection phenomenol- 
ogy of composite models arise because the dark matter 
has a finite size Aj 1 3> m^T. 1 . The cross section is sup- 
pressed by an effective form factor when a neutral bound 
state interacts with momentum \q \ <C Ad |26j . States 
with nonzero spin have multipole interactions with the 
field. These moments vanish for states with zero spin; 
scalar states that can only couple through the charge- 
radius and polarizability interactions are the dominant 
scattering mechanisms. For the dark pion scattering off 
the SM, the charge-radius interaction dominates over the 
polarizability interaction, which is suppressed by an ad- 
ditional factor of the mixing parameter e 2 . 

The charge-radius is the effective size of the tt^ 
probed by the dark photon. In the limit of small mo- 
mentum transfer \q\ -C Ad, the wavelength of the dark 
photon is too long to probe the charged constituents of 
the composite state and the scattering rate is suppressed. 
Elastic charge-radius scattering cannot be the sole con- 
tributor to the direct detection signal due to constraints 
from current null experiments. However see |241 125] for 
examples on how form factors can reconcile DAMA with 
the null experiments. 

The dominant scattering is inelastic and there is a 
subdominant elastic component that accounts for a frac- 
tion of the total scattering rate. Specifically, the differ- 
ential scattering cross section is 



da 

dEn. 
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m N E R 
2m^ A 5m 
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where is the mass of a nucleus with charge Z recoil- 
ing with energy E R , c? n ,c"j = 1 of Sec. II C| and 



dap 8Z 2 arriN 1 



dEr 
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feff (l + 2m N E R /m 2 



(18) 



1 This only applies to the field strengths, , not the gauge 
potentials, A*f, whose interactions are constrained by gauge in- 
variance. 



The scattering operators couple the dark matter states 
coherently to the nuclear charge, and the Helm form fac- 
tor accounts for loss of the coherence at large recoil 



|^Holm(^fl,)| 2 — 



3ji(kko) 



(19) 



where s = 1 fm, tq = yV 2 — 5s 2 , and r = 1.2/1 1 / 3 fm 

m 

The differential cross section depends on the con- 
finement scale Ad = \J m^ A 5m/ k, the mass of the dark 
photon m,A d , and the couplings of the effective theory 



/off : 



(20) 



where k is the 0(l/N c ) constant defining the mass dif- 
ference from Eq. |4j 

To determine the preferred region of parameter 
space for CiDM models, a global x 2 analysis was per- 
formed that included the results from all current direc- 
tion detection experiments. This procedure is outlined 
in [3SJ|23] and is summarized here. The differential scat- 
tering rate per unit detector mass is 



dR 
dEn 



-^^[d*vf(v + v e )v^, (21) 
m^ d m N J dE R 



where po = 0.3 GeV/cm 3 is the local dark matter den- 
sity and v e is the velocity of the Earth in the galactic 
rest frame. There are significant uncertainties in the 
dark matter velocity distribution function f(v), and con- 
straints on the particle physics model can vary wildly 
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FIG. 1: 95% contours in m 7rd — Ad parameter space for df = 
0.00,0.06,0.07,0.08. For this figure, k = 1/4 and m. 4d = 1 
GeV. The dashed lines show contours of Sm in keV. The inset 
shows the 68% confidence regions for df = 0.00, 0.04, 0.06 
for the same k and m,A A . The colors correspond to 8j> = 
0.00 (blue), 0.04 (teal), 0.06 (magenta), 0.07 (yellow), 0.08 
(green). 
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FIG. 2: 95% confidence limit regions of m 7Td — f e g for Qf = 
0.00,0.06,0.07,0.08. For this figure, k = 0.25 and m Ad = 
1 GeV. The inset shows the 68% confidence regions for 9f = 
0.00, 0.04, 0.06 for the same k and rriA A - 



depending on the particular choice of benchmark halo 
model. To find the full scope of allowed CiDM models, 
we marginalize over a parameterized velocity distribu- 
tion function of the form: 



f(v) oc exp 



exp 



(22) 



where the parameters are constrained to be within 

200 km/s <v < 300 km/s 
500 km/s < v csc < 600 km/s 

0.8 < a < 1.25. (23) 

These values are motivated by observational constraints 
[301 131] and analytic approximations to the Via Lactea 
results 

The global \ 2 fit is performed by marginalizing over 
the six unknown parameters of the dark matter and halo 
model: m Vd , 8m, f e s, Vq, v esc , and a. The measurements 
used in the \ 2 fit are the first twelve bins of DAMA's 
modulation amplitude, as well as a single high energy 
bin from 8 keVee to 12 keVee [U [35]. In addition 
to DAMA's signal, the dark matter predictions are re- 
quired to not supersaturate any observation from null 
experiments at the 95% confidence level. The null exper- 
iments included in the analysis are: CDMS [fS"l ITB], H7]. 
ZEPLINII [35], ZEPLINIII [37], CRESSTII [35J [3J], 
and the new XENON10 inelastic dark matter analysis 

001 HI]. 

The lcr and 2er allowed regions in the m ld — Ad and 
m TT d ~ fcS spaces are shown in Fig. [T] and [2] respectively. 
The minimal \ 2 nas a value of 4.61 and the correspond- 
ing point is listed as model CiDMi in Table [T] The 1 a 
and 2 a regions are set by 



where Ax 2 (ler) = 7.0 and Ax 2 (2ct) = 12.6. There- 
fore, the 68% and 95% regions are set by requiring that 
X 2 < 11.6,17.2, respectively. As the fraction of form 
factor elastic scattering increases relative to the inelastic 
contribution, the allowed regions in Fig.[T]and[2]each sep- 
arate into two. At 95% confidence, dark matter masses 
with m Vd > 200 GeV correspond to "slow" velocity dis- 
tribution functions where vq < 225km/s and a > 1.15. 
One benchmark model is shown in Table [T] as CiDM3. 
However, the correlations between the dark matter mass 
and the velocity distribution parameters are far weaker 
in the low mass region (m„ d < 200 GeV). 

The mass of the dark photon is related to the mixing 
parameter e as 



2V2 



9df 2 B 



VhI 2 



(25) 
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For theories with dominant inelastic scattering, / e ff ~ 
O(700 GeV) and m OTd - O(100 GeV) to satisfy both the 
DAM A and null experiments. Therefore, keeping ~ 1 



O(10~ 4 ) 



1 GeV 



(26) 



which corresponds well with the results of the \ 2 global 
fit. Fig. [3] shows the lcr and 2cr regions in the rriA d — e 
parameter space allowed by all current direct detection 
experiments. A benchmark value of yu = 1 is chosen; 
the contours shift to larger e for smaller Yukawa cou- 
pling. 

Light Ad alter the fit to the DAMA spectrum be- 
cause the propagator suppresses high momentum scat- 
tering events. The momentum transfer needed to explain 
the highest energy bin with a statistically signifiant an- 
nual modulation rate in the DAMA spectrum [Er = 5 
keVee) is 



i2m I E R 
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120 MeV. 



(27) 
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(24) 



If the mass of the dark photon is less than 120 MeV, its 
propagator in (18 1 suppresses the scattering rate in the 
high energy bins. The suppression of the high momen- 
tum transfer events can be compensated if the mass split- 
ting, 5m, grows larger; however this eventually forces 
f~fi to become large, increasing the allowed values of e. 
These effects are shown in Fig. [3] A low bench- 
mark model is shown as CiDM4 in Table. [I] The follow- 
ing section presents constraints on the allowed parameter 
region arising from indirect and direct searches for the 
dark photon. 



IV. SEARCHES FOR THE DARK PHOTON 

The dark photon communicates with the Standard 
Model through kinetic mixing and experimental bounds 



5 















72 GeV 


75 GeV 


234 GeV 


162 GeV 


5m 


109 keV 


105 keV 


91 keV 


126 keV 


A 

A d 


Y( ( Me V 


1 77 A ft ^\ T 

lit Me V 


zyz MeV 


zoo MeV 


/eft 


70 n ~\ 7 

too LreV 


o4o OeV 


OOO CjreV 


zoo Lrev 


£ 


O. i X 1U 


^ n v 1 n — 4 


9 i v i n — 3 

Z . _1_ X 1U 


O.O X 1U 


vriA , 


1 GeV 


1 GeV 


1 GeV 


60 MeV 


6j> 


0.00 


0.04 


0.04 


0.06 


V() 


272 km/s 


273 km/s 


202 km/s 


280 km/s 


Vesc 


510 km/s 


501 km/s 


558 km/s 


501 km/s 


a 


0.86 


0.82 


1.30 


0.98 


x 2 


4.6 


6.2 


12.8 


9.9 



TABLE I: Four benchmark models showing different regions 
of parameter space. CiDMi corresponds to the best-fit point. 
QDM2 shows a representative mixture of inelastic and sub- 
dominant elastic scattering. GDM3 shows the larger mass 
window with slow halo parameters. C1DM4 shows a light 
rriA A model. 



on these interactions arise from tests of QED. The most 
model-independent bound comes from the virtual ex- 
change of the A d between SM fields. The best limits arise 
from the constraints on the magnetic dipole moments of 
the fj, and e [42]. The constraints can be expressed as 



ysis is beyond the scope of this paper and will be per- 
formed in [43] . This article uses a bound on e of [44] [45] 



e<lx 1CT 2 . 



(31) 



This constraint becomes more important than the muon 
g - 2 limit when m Ad > 250 MeV. 

The width of the Z° is altered by the presence of 
the dark sector. The interaction between the Z° and the 
dark current in the canonically normalized mass eigen- 
state basis is 
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to lowest order in e and m 2 Ad /rn 2 z0 - The width of the Z° 
decay into the dark sector is 
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Any additional Z° decay mode cannot have a branching 
ratio of more than 0.18% [46]. This sets a limit on the 
dark photon mass of 
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which is never a constraint. 
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where 
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(29) 



Fig. [3]shows that g — 2 limits are most important at low 
m Ad and large e. 

For larger values of m Ad , constraints arise from 
precision electroweak interactions, which depend on the 
gauge terms of the Lagrangian. The kinetically-mixed 
C/(l)d only alters the neutral currents and the La- 
grangian for this sector is 
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The precision electroweak constraints have not been per- 
formed for dark photons with masses between 1 GeV and 
100 GeV. In addition to oblique corrections, there is 
a non-oblique correction coming from the contribution 
of the dark photon to precision electromagnetic observ- 
ables, such as differential Bhabha scattering. A full anal- 



The allowed parameter space shown in Fig. [3] can be 
further constrained by searches for the direct production 
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FIG. 3: The 95% limits on m. 4d - e with y H = 1 for 6 r = 
0.00 (blue), 0.06 (magenta), 0.08 (green). The dark gray 
regions are excluded by limits on the g — 2 of the electron 
and muon (top left) and fixed-target experiments (light vtia a 
and moderate e). The light gray region shows limits from the 
BABAR T(3S)— > jh + search; the direct search limits are 
model-dependent and must be interpreted on a case-by-case 
basis. 
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of the A d @7J US] HO |5D1 EE]. In this section, we out- 
line the prospects for such searches and the challenges of 
translating the experimental bounds to theoretical con- 
straints in composite models. 

If the dark photon is the lightest state in the dark 
sector (mA d < Ad), then it will decay directly to the 
SM. Such a light Ad will be dominantly produced from 
e + e~ — > j Ad and from the decays of the dark hadrons. 
Once produced, the dark photon will decay promptly; 
for example, 
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(35) 



for the benchmark model CiDMi. Hadronic decays are 
also allowed, but are subdominant to the lepton decays, 
except near resonances [4"T] . 

When the dark photon is heavier than Ad, it can 
either decay to dark mesons or directly to SM leptons. 
However, the coupling of the dark photon to the electro- 
magnetic current is suppressed by a factor of e relative 
to the coupling to the dark quarks 



C 



dd 



c 



dcm 



(36) 



In this limit, the branching fraction into SM leptons is 
negligible: 
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for N c = 4 and the benchmark model CiDMi. The dark 
photon preferentially decays to the dark quarks, which 
first parton shower, then hadronize, and finally cascade 
decay back to SM particles. 

The most common mesons formed in the hadroniza- 
tion process will be the lightest in the spectrum. There 
are no light pseudo Goldstone bosons in this theory be- 
cause there is only one light quark. Therefore, the light- 
est meson is rj' d (0 h ), in analogy with QCD. Using the 
SM rf as a prototype and the SM ao meson as a typical 
hadronic state, the mass of the r]' d is estimated to be 
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The f]' d becomes light in the large N c limit. For N c > 10, 
the dark photon can decay to the r\' d . However, rj d 
is cosmologically stable because it has a chirality sup- 
pressed decay to electrons and primarily decays via a 
loop-induced process to two leptons, which dominates 
over the four body decay rj' d — ► A d A* d |47l|52]. The decay 
width is suppressed by an additional factor of (Ad//^) 2 
relative to the <f> decay mode of [171 1521 [53] because the 
r]' d has to mix with the <fi to mediate the decay. The 
resulting decay width is 



T( v ' d - e+e") 



e 4 a 2 m 2 m r/ Aj, 1 

C 10iVars' ( j 



The cosmological relic abundance of r)' d is sufficiently 
small to make up a small fraction of the matter density 
of the Universe. 

Unlike the rf d , the next-lightest meson, tj d 
can have prompt decays. The mass of the cod also be- 
comes small in the large N c limit. The uid will decay to 
SM leptons by mixing with the dark photon. Approxi- 
mating the mixing angle by 



rrr. 



the decay width is 



r(wd 



■o 



e 2 ag 2 d c 2 g < 
3 m\ 



ac 2 g A d 



3 / e 4 ff 20 m 



(40) 



(41) 



for the CiDMi benchmark point in Table [T] Therefore, if 
the u>d is produced, it will decay to two leptons with a 
long displaced vertex. 

The dark photon will decay promptly to leptons 
if m Ad < A d and both the BABAR [53] and CLEO 
[55] searches for T(3S)^ 7^ + /i - may be used to set 
bounds. However, when the dark photon decay chan- 
nels are closed, the muon decay channels are often closed 
as well because Ad ^ 2m M . The estimated bounds from 
these searches are shown in Fig. [3] It may also be pos- 
sible to use the T(IS') — * 7 + X, where X is invisible, 
when the A^ decays outside the detector [5B] . 

The best chance of discovering the dark sector is 
by directly producing the ^4d at low-energy lepton col- 
liders. BABAR has recently searched for the A^ in the 
Ai channel |57j and future work is being pursued at a 
myriad of experiments [581 159] , The searches are com- 
plicated because the decay of the A d back to leptons is 
suppressed by the factor in Eq. [37] To gain efficacy, it 
is necessary to use the decay into the dark quarks. For 
m^/Ad ^ O(10), the number of dark hadrons is moder- 
ate. Because the fluctuations in the number of hadrons is 
non-Gaussian, the cost of fragmenting to two w d mesons 
is not limiting and it is possible to set limits using the 
BABAR and CLEO searches. When m^/Ad > 1, there 
can be a large number of dark hadrons in the decay prod- 
ucts of the Ad, and an inclusive, multi-lepton search is 
necessary. When the dark photon decays through the 
hadronic channel, the analysis becomes more challeng- 
ing because it is necessary to know how the dark par- 
tons fragment into dark hadrons and then decay down 
to the Standard Model. Setting limits on this model is 
beyond the scope of this work because of the significant 
uncertainties in the hadronic spectrum. 

In addition to low energy searches for the decay 
products for the Ad, high energy colliders provide a use- 
ful laboratory. LEP-I can search for rare decays down 



to the Br(Z°) < 10~ 5 . From Eq. |33l this corresponds 
to masses of the Ad ^ 4GeV. When the Ad decays with 
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a mass mA d ^ Ad, it decays into a pair of dark quarks, 
^l^l, and proceeds to shower and hadronize. Future 
studies of LEP2 are needed to determine the relevant fi- 
nal states and the procedure necessary to set limits on 
the hetrogeneous final states. 



V. DISCUSSION 

This article introduced a composite inelastic dark 
matter model with dominant inelastic scattering off of 
nuclei and a subdominant elastic scattering component. 
The subdominant elastic component is a signature of the 
symmetry structure of the model and is a critical feature 
to measure. It was found that parity violating effects can 
be nearly maximal, with 



e, 



a7 



<0.' 



(42) 



Discovering the elastic subcomponent will place a lower 
limit on tol and could sharpen the Standard Model's 
strong CP problem and flavor structure. 

Directional detection experiments, which measure 
both the energy and direction of recoiling nuclei in de- 
tectors [BHl ED E2] , can be used to distinguish the elastic 
and inelastic scattering components by looking for large- 
angle scattering events [25] . There is an upper bound on 
the allowed scattering angle, which depends on the types 
of interactions that are allowed. In particular, 



COS7 n 



(43) 



where 7 is the angle between the direction of the Earth's 
velocity and the recoiling nucleus in the lab frame and 
Vmin is the minimum velocity to scatter at a given recoil 
energy, 




elastic 
inelastic 



,(44) 



and n is the reduced mass of the dark matter-nucleus 
system. Inelastic scattering events have a much smaller 
cos7 max than elastic scattering events and these two 
types of interactions can be distinguished with the next 
generation of directional detection experiments j63j EH 

E5JH3E7]. 

The best fit for CiDM had Sm ~ 100 keV and m n<i ~ 
70GeV. This leads to an estimate of 



A, 



150 MeV 



(45) 



for the dark sector confining scale. The lower 
bound is Ad <^ 70 MeV / y/n and the upper bound is 
A d < 240MeV/V^ (for Fig. [l] a value of k = 0.25 
is used). This indicates that the dark matter form 



factor may be important in shaping the higher energy 
bins of DAMA. Applying a dark matter form factor 
might therefore change the allowed parameter range. For 
Coulombically-bound dark matter, the form factor can 
be found by Fourier transforming the hydrogenic wave 
functions to get 



1 



1 + 9 2 ?DM ' 



(46) 



where 



DM 



= ^d m L is the Bohr radius of the 



Coulombically-bound state and Ad is the 't Hooft cou- 
pling evaluated at the Bohr radius. The strongly inter- 
acting form factors can be estimated by extrapolating 



DM 



Ad and behave similarly to having rriA d 



DM- 



CiDM models have sub-components to the dark 
matter that are not the tt^. There are roughly three 
classes of particles: the /?d, multiple heavy quark mesons 
(e.g., ^ ' H ty H ^> L ^> L states), and baryons (e.g., ^ h 1 • • * h 
states) . The relative populations of these states is deter- 
mined by the interactions in the early Universe [5j [68] . 
Detecting the latter two classes of particles would be a 
clear indication of composite inelastic dark matter. The 
signature will be striking because the mass of the dark 
matter subcomponents would be near-integer multiples 
of the 7Td mass, ending at N c m^ d . 

The collider signatures of CiDM are challenging be- 
cause many of the decay products have long lifetimes 
and give rise to extremely displaced vertices or missing 
energy. To interpret the results from current e + e~ col- 
liders, it is necessary to have better estimates for the 
dark hadron multiplicity distributions from dark pho- 
ton decays. If there are dynamics that stabilize the vev 
of both the Standard Model Higgs and the dark sector 
Higgs, then it is possible to produce dark sector states at 
the Tevatron and LHC in the decays of electroweak scale 
particles. The phenomenology of these events will be 
similar to that of hidden valley theories, with the major- 
ity of parameter space giving rise to extremely displaced 
vertices |691 170L 171] . Beam dump experiments are ide- 
ally suited for identifying leptons from particles with a 
finite lifetime and provide the best prospect for discover- 
ing the dark sector through direct production. Proposals 
for these experiments are presently underway |59j . 
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